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Abstract 

We calculate and compare the response of light wave interferometers and matter 
wave interferometers to gravitational waves. We find that metric matter wave in- 
terferometers will not challenge kilometric light wave interferometers such as Virgo 
or LIGO, but could be a good candidate for the detection of very low frequency 
gravitational waves. 
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1 Introduction 



The extreme sensitivities required to detect gravitational waves still represent 
a challenge. LIGO interferometers 1 have reached their nominal sensitivity and 
have begun scientific runs giving upper limits of gravitational waves in various 
cases : gravitational bursts, stochastic background and periodic waves [1,2,3]. 
Virgo 2 is now in the commissioning phase; its sensitivity should be better 
than LIGO's at low frequencies (around 10 Hz). 

However despite important technological improvements in many domains, the 
present detectors will not be able to observe the very low frequency sources 
of astrophysical interest. The LISA project [4], a space-based laser interfer- 
ometer with 5 millions kilometers arms, is presently the best known solution 
to improve the detectors capabilities at very low frequencies. The challenge is 
impressive, we thus believe that new technologies will be necessary for further 
developments in the future. 

In 2004, Chiao and Speliotopoulos [5] proposed to use matter-wave interfer- 
ometers as gravitational wave detectors (the so-called MIGO, Matter- wave 
Interferometric Gravitational- wave Observatory). They claimed that MIGO 
could reach the same sensitivity as LISA or LIGO in a much more compact 
configuration. Their results are discussed in papers from Roura et al. [6] and 
Foffa et al. [7]. The various authors find divergent results. These differences 
originate from different physical interpretations of the various coordinate sys- 
tems which were used and from the different boundary conditions which were 
assumed. 

1 http://www.ligo.caltech.edu/ 

2 http://www.ego-gw.it/ 
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In the present paper we consider various matter- wave interferometers. For each 
of them, we estimate the order of magnitude of the phase shift due to a periodic 
gravitational wave. We compare the sensitivity of matter wave interferometers 
(MWI) and light wave interferometers (LWI). Several boundary conditions 
are considered : fixed or free mirrors. In each case we clearly point forward 
the corresponding comoving coordinates. We emphasize the role of the three 
mostly relevant parameters, i.e. the atomic flux, the atomic speed and the 
number of atomic bounces in the arms of the interferometers when similar to 
Fabry-Perot cavities in optics. 

In the present paper, we do not consider the technical feasibility. Our goal 
is just to compare MWIs and LWIs in order to become familiar with MWIs 
and to discover the conditions necessary for matter wave interferometers to 
be useful. Once these conditions are put forward it might be possible to have 
an opinion on the future feasibility of such detectors. For this next step, an 
expertise in high-tech experiments, cold atoms and matter wave interferometry 
is crucial, as well as a high dose of optimism. 

Now, before we start, let us recall that a gravitational wave is a perturbation, 
h^ v , in the local metric : 



where r]^ is the usual Minkowski metric and the greek indices run from to 
3. In the sequel h refers to the order of magnitude of the biggest j/i^j (i.e. 




KA ^ h « !)■ 
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2 The quantum phase 



In this section we consider the matter wave associated to a particle of mass 
m (for a light wave m = 0). We investigate the perturbation of the quantum 
phase due to the perturbation h^ u of the metric. 

First we consider the propagation of the wave between two mirrors or two 
beam splitters A and B. The quantum phase is developed up to the first 
order relatively to ft : (ft = (p + 5(f) where (fi is the unperturbed quantum phase 
while 5(f> is the perturbation of order h. The phase is assumed to be a solution 
of the eikonal equation : g^ v <i)^(t>,u — {mc/h) 2 , where ft is Planck's constant 
and c the light velocity in vacuum. Therefore the unperturbed solution is 
4> = k^ + est with = = rj^p" /ft, where p u is the unperturbed 

particle momentum. For a matter wave p v = r ymc(l,v/c), where v is the 
group velocity of the wave and 7 = . Let us notice that the atom 

energy is E = ftck . 




Fig. 1. The matter wave propagates between the two optical elements A and B 
(dashed line). The unperturbed trajectory is the straight line. 
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The wordlines of the mirrors or the beam splitters A and B are perturbed by 
the GW. Let call x\ = x^(t) and 

B — ^ B 

(t) their trajectories. Knowing 
and t B (the arrival time of the atom, we can deduce t A (its departure time) : 
x b ~ x a = a ^-> where a is a constant. This is illustrated on Fig. 1. We define 
[0]^ = (f> [ig (t B )] — <p [x A In order to calculate [0]^ we use the method 
described in [8] and [9]. The solution reads : 

K = [Ml + Ml (2) 

where : 

['PoIa = kn x B ~ kn x A (3) 

W]* = ^£V«"§ (4) 

The integration in (4) is performed along the unperturbed trajectory (the 
straight line of Fig. 1). 

The phase is a scalar, therefore [4>} A is independent of the coordinate system, 
although the decomposition (2) is not. A first order coordinate transformation 
results in a first order change of the functions x^ A B {t) and a first order change 
of h^ u in (1). Both effects have to be into account. Then, contrary to the 
claim of Chiao & Speliotopoulos in [5], one can show explicitely that the 
result is invariant [6]. In the sequel, for the sake of simplicity, we will choose 
a coordinate system where the optical elements are at rest, i.e. xa and x B are 
independant of the time coordinate x°. 
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3 The free interferometer 



3.1 The Einstein Coordinates (EC) 



In the weak-field approximation, in empty space, the Einstein equations are 



RfA,u 2^ {h<TT,tiv ~\~ h/j,ii,aT h uu ^ T h^ Ti0 - u ) (5) 



In harmonic coordinates the condition d u h^ = \r)pd u h° holds true. The Ein- 
stein equations (5) read : 



□V = ( 6 ) 

Gravitational waves (GW) are solutions of such a propagation equation. Far 
away from the sources, plane waves are solutions of (6). Some more constraints 
on h^y (i.e. <9j/i* = and h\ = 0, where latin indices run from 1 to 3) define 
an unique coordinate system : we will call it the Einstein Coordinates (EC). 
To simplify the problem, we assume that the GW is propagating along the 
z = x 3 axis. The only non zero components of are h rs = h rs (z — ct), with 
r, s = 1 or 2, then : 



ds 2 = T]^ dx 11 dx u + h rs dx r dx s (7) 



In the slow motion limit, the space coordinates x l of a point mass particle 
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satisfy the equations 



(iV -h)^f + 0(h 2 ,v 2 h) with 



dt 2 



dx l 



dt 



< c 



where () = —(). 



It is easy to see that a point mass initially at rest remains at rest in the EC 
during the passage of a GW. 



3.2 The Michelson-Morley interferometer 

We study a Michelson-Morley configuration with the arms along the x 1 and 
The optical elements are supposed to be free of any constraints (free 
falling in LISA or fixed on super attenuators in Virgo). Their spatial EC 
remain constant with time. 



L 






L 




<- -> 


k 


^ A^> = <j) 2 ~ ^ 



Fig. 2. Michelson interferometer. 



We assume that h rs is a sinusoidal function where is the angular frequency of 
the gravitational wave. Moreover we assume that the size of the experimental 
device along the z axis is much smaller than the gravitational wavelength A = 
2ttc/Q. These conditions results in hu (z — ct) = —h 2 2 (z — ct) = h + ■ sin (Qt) 
and hi2 {z — ct) = /121 (z — ct) = h x ■ sin (Qt + (p x ) where h + , h x and (p x are 
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three constants. Then, using the notations of Fig. 2 and eqs. (2) and (7), the 
phase difference between the two arms of the interferometer is : 




(9) 



where V = c for a light wave and V — vq, the initial group velocity, for a 
matter wave. A is the wavelength; for a matter wave, it is the De Broglie 
wavelength : A = 2nh/m^VQ. 

When VtL/V < 2ix the amplitude of A0 in (9) is : 



This formula is well known for LWIs such as Virgo or LISA [10, p. 54]. It holds 
true for a MWI. However, for a LWI, the condition QL/V <C 2ir corresponds 
to L < A while it implies L < — • A for a MWI. 



3.3 Light Wave Interferometers versus Matter Wave Interferometers 

In Virgo or LISA, the laser source is a Nd:YAG infrared laser with A ~ 1CT 6 m. 
Decreasing the laser wavelength in order to increase the amplitude A(f> of the 
phase difference in (10) raises technical difficulties (laser stability, mirrors 
quality). On the other hand it seems easier to achieve in a MWI : for instance 
in [11], the wavelength is 16 pm. As a consequence, the phase difference in a 
MWI can be much larger than in a LWI, and the sensitivity appears to be 
much better. Unfortunately, this conclusion does not hold anymore when one 
considers the fundamental limit. 



A0 = Aix\h 




(10) 



c 
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The signal is ultimately limited by the shot noise. The minimum phase differ- 
ence that can be detected is : 



A0 = (11) 

2VNt 

where N is the particle flux and t the integration time. 

The atom flux in the experiment described in [11] is only 70 s _1 which is very 
low but in reference [12] the atom flux is 10 11 s _1 for a velocity of 290 m.s -1 (i.e. 
A ~ 10 pm). Such a flux remains however very small compared to the photon 
flux in Virgo which is of order 10 23 s" 1 [13]. Therefore one can show that 
metric MWIs have to be considerably improved to compete with kilometric 
LWIs. 

However, the LISA effective flux can be estimated of order 10 s s" 1 [4, p. 60]. 
It is much lower than the Virgo effective flux, and even lower than the better 
MWI fluxes. Therefore, it seems easier to achieve with a MWI the sensitivity 
of LISA than the sensitivity of Virgo. 

Assuming that the detection is only limited by the shot noise and that the inte- 
gration times are the same, one obtains the conditions for similar sensitivities 
of MWIs and LWIs from formulas (10) and (11) : 



2nhL 



lw 



lw 



[12) 



where the subscripts mw and i w denote respectively the characteristics of the 
MWI and the LWI. 
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Fig. 3. Required caracteristics of a Fig. 4. Required caracteristics of a 
MWI necessary to reach the sensitivity MWI necessary to reach the sensitivity 
of Virgo. of LISA. 

On Fig. 3 and 4 we represent the required characteristics 3 of a MWI necessary 
to reach the sensitivity of Virgo (Fig. 3) and that of LISA (Fig. 4). The 
characteristic point of the MWI described in [12] has been represented on 
both figures. 

It appears clearly that the atom velocity has to be very high in order to 
reach Virgo sensitivity. However it would be a difficult challenge to keep the 
coherence and to separate the matter ray at such velocities. An idea could be 
to accelerate atoms inside the interferometer arms after separating the matter 
ray, and decelerate it before the reflections. For example, in GANIL 4 , one can 
obtain an ion ray at half the light velocity with fluxes up to 2,6. 10 13 s" 1 . In 

These curves are drawn from formula (12) for the caesium mass. vq : Zy mw and 
N mw are respectively the initial atom wave group velocity, the MWI arm lenght 
and the atom flux. The characteristic point corresponds to the MWI described in 

[12]. 

Grand Accelerateur National d'lons Lourds 
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order to decrease the velocity a possibility could be to increase the matter flux 
or the mass of the atoms [14] or to increase the arm length. However a good 
compromise cannot be achieved nowadays. 

On the other hand, it seems easier to reach LISA sensitivity since the required 
atom velocity is much lower. The characteristic point of the MWI described 
in [12] corresponds to a kilometric interferometer. A one meter MWI with 
~ 1, 000 round-trips in each arm would play the same role. It would be similar 
to the Fabry-Perot cavities in Virgo like interferometers. 

4 The rigid interferometer 

4-1 The Fermi Coordinates (FC) 

In usual matter wave interferometers the mirrors and the beam splitters are 
"rigidly" bounded [11,12,15]. The concept of rigid body is not a relativistic 
one, however if the experimental set-up is much smaller than the gravitational 
wavelength (L C A), one can introduce Fermi like coordinates (FC) comoving 
with the matter of the set-up supposed to be "rigid" [10, pp. 39-47]. We will 
assume here that a rigid body has fixed coordinates in the FC. 

Lasers are used as mirrors and beam splitters in the experiments [12,15] : it 
does not add a major perturbation to the phase difference. Indeed, in FC, the 
description of the physical phenomena (including the Maxwell equations) is 
very similar to the special relativistic one [10, p. 52]. 

We choose the center of mass of the interferometer as the origin of the Fermi 
reference frame. The whole device is supposed to be free of constraints, so that 
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the worldline of its center of mass is a geodesic. The metric can be derived 
from the general formula of Manasse & Misner [16]. It is convenient here to 
derive the coordinate transformation from EC, x a , to FC, X a , in order to 
link the movement of the optical elements in the two coordinate systems. This 
coordinate transformation has been derived with a general method that will 
be discussed in detail in a future paper. In order to find the metric up to the 
second order relatively to the X\ one needs the coordinate transformation up 
to the third order : 

x r = X r _ l h r X s _ l^'Z 2 + 0(£ 4 , k 2 ) 

(13) 

x a = X a - \h rs X r X s - l'h rs X r X s Z + 0(£ 4 , h 2 ) 



where r, s = 1 or 2, a = or 3, Z = X 3 , £ ~ sup \X % \ and c = 1 (geometrical 
units). Then, the metric is : 



ds z = r]^ dX^ dX u + -h rs X r X s [dT z + - dZ 2 - -h rs X r X s dT dZ 



2 V 



+-h rs Z dX r (2X S dT -X s dZ + \zdX s ^j + 0(£ 3 , h 2 ) 



(14) 



where T = X°. This metric is in agreement with reference [17] where Fortini 
& Gualdi obtain the metric up to any order, although they do not derive the 
coordinate transformation that we consider here. 
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In the neighborhood of the plane Z — 0, one obtains 



x r = x r - \h r s X s + 0{£ 4 ,h 2 ) 

(15) 

x a = X a_ lh rs X r X S + 0(£ 4 , h 2 ) 

and : 

ds 2 = r]^ dX» dX v + \h rs X r X s dT 2 + 0(f, h 2 ) (16) 

The coordinate transformation (15) was first found by Ashby & Dreitlein in 
reference [18]. However, contrary to the claim of the authors, the metric (16) 
(that can be derived from (15)) is a Fermi metric in the sense of Manasse & 
Misner [16] only in the plane Z = 0, and not in the whole space. 

Now that we have the metric tensor in the FC, we will compute the phase 
difference for several rigid interferometer. 

4-2 The Michelson-Morley configuration 

First we study the Michelson-Morley configuration of Fig. 2, but now we as- 
sume that the coordinates of the optical elements are constant in the FC. 
From eqs. (2) and (16) we find : 

A0 = Anh + ■ y • fl - 7Tr sin T7"V sin ^ ( 17 ) 
A V i ILi V J 

with the notations of section 3.2. This formula is in agreement with the one 
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obtained in [6]. The assumption L <C A implies VtL/c <C 2ir. Therefore, for a 
LWI, formula (17) reduces to : 

A0 lw ~ (18) 

This result is well-known : if the arms of a LWI were rigid there would be 
no signal. The situation is different for a MWI where we can consider two 
different regimes [6] : 

• L < — • A (i.e. QL/V < 2ir) : formula (17) reduces to A0 mw ~ 0. This 

c 

regime occurs when the flight time of an atom in the interferometer is much 
less than the GW period. In this case there is no signal. 

• — -A<I (i.e. QL/V 3> 2ir) : one can assume V/VlL <C 1. This regime 

c 

occurs when the flight time is much longer than the period of the GW. In 
this case the amplitude of the phase difference in (17) reduces to : 

A0 mw ~ 4nh + ■ j (19) 

This new specific regime has no equivalent with free interferometers nor 
with a rigid LWI. One can notice that the present amplitude given by (19) is 
similar to formula (10). Unfortunately in the present regime the sensitivity 
is limited by the condition t> flL. 

4-3 The Ramsey-Borde configuration 

We study now a MWI with a Ramsey-Borde configuration [11,12,15] repre- 
sented on Fig. 5. We define X = X 1 and Y = X 2 . Atoms move along the 
X axis and we assume that Vq <^ c. At the output of the interferometer, we 
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atoms- 




Fig. 5. Ramsey-Borde interferometer in the X-Y plane, 
calculate from eqs. (2) and (16) the phase difference between the two beams : 



A0 = A0 + + A0> 



(20) 



with : 



2 

-4irh + — sin ^ tan 9 
A 



sin ip 

cos (fit + ip) H — cos (fit) 



f21} 



A0 X = — 47r/i x — cos ■i/' tan ^ 
A 



sin (fit + ip x — ip) 



sin?/; 



tan ip cos (fit + <^> 



(22) 



QL 

where 9 is the separation angle and ib = . 

2v 

The lasers give to the atoms a velocity Av = vq tan 9 along the Y axis. In 
order to achieve a high value of 9, one must communicate to the atom a high 
momentum during the interactions with the photons. This is very difficult 
because this has to be done without any loss of coherence. The momentum 
of the photons being hk, one can imagine to transfer mAv = N x hk to the 
atom. Significant results, N ~ 1000, have already been obtained along this 
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line of research [19,20]; they still remain insufficient today. 

We assume now that 9 ~ 7r/4 so that tan 9 ~ I. Here again two regimes can 
be considered : 

• L <C — ■ A (ie) ib = <C 2ir : in this case 

c 2v 

Acb ~ -4nh x ■ ^ • sin (fit + ip x ) (23) 
A 

This expression results in an amplitude already given, formula (10) above. 
Therefore, Figs. 3 and 4 can be used for the discussion of the experimental 
design of the interferometer (with h + — > h x ). 

• — -A<L (ie) 2n <^ = ib : expression (20) gives 
c 2vq 



A0 ~ —An ■ y • [h+ sin ib cos (Qt + ip) + h x cos ib sin (fit + {p x — ib)] (24) 
A 

The phase difference A0 is a periodic function of the time whose ampli- 
tude displays the same order of magnitude that was previously put forward 
in expressions (10) and (23). Here, however, one can choose the value of ib 
in order to measure h + or h x . This is a positive point but the sensitivity 
remains however limited by the low values of vq in this regime. 



5 Conclusion 

In this paper we considered matter wave and light wave interferometers de- 
signed to detect gravitational waves. This has been considered already in the 
literature in some special cases [5,6,7] where several different claims corre- 
spond to different assumptions (more or less justified) about the coordinate 
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systems comoving with the interferometer. 

We considered "free" and "rigid" interferometers. The comoving coordinates 
have been chosen as Einstein coordinates in the first case; it seems that this 
choice raises no discussion. The comoving coordinates have been chosen as 
Fermi-like coordinates in the second case. This choice is very natural because 
one can show that the first approximation of the equations of a continuous 
medium are ordinary non relativistic equations with just an extra gravitational 
force density (the mechanical detectors of gravitational waves are precisely 
based on such equations). 

We have considered only plane interferometers orthogonal to the propagation 
of the gravitational waves. If we change the orientation of the interferom- 
eter relatively to the gravitational wave, the sensitivity is slightly modified 
but the orders of magnitude remain the same. Moreover we only considered 
periodic waves and no pulses of gravitational waves. We believe that the com- 
parison between light-wave interferometers and matter-wave interferometers 
is not deeply affected by such a simplification. 

In the cases that we studied we obtained an estimation of the sensitivity of 
matter- wave interferometers. We especially considered the shot noise limit, 
however, it is the thermal noise which most troublesome at the present mo- 
ment. The answer of this problem could be the construction of a compact (one 
meter) interferometer which could be cooled at very low temperature. In order 
to estimate roughly the thermal noise we consider that the interferometer is 
fixed on a bench of mass M ~ 500 kg which displays an eigenfrequency of 
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order u ~ 10 4 s 1 . Therefore, following reference [13], one finds the limit : 



where ks is the Boltzmann constant. With a quality factor Q ~ 10 , and a 
temperature T ~ 10 -2 K one finds h min ~ 10~ 20 for Q ~ 10~ 2 s _1 . 

Now if a clear conclusion had to be taken from the previous estimations, we 
would claim that in the future, compact, very low temperature matter- wave 
interferometers will not be a serious challenger to high frequency detectors of 
gravitational waves (such as Virgo or LIGO) but to LISA. Of course major 
improvements still remain necessary today. 
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